Second sound spectroscopy of a nonequilibrium superfluid-normal interface 
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An experiment is proposed to test a previously developed theory of the hydrodynamics of a 
nonequilibrium heat current-induced superfluid-normal interface. It is shown that the interfacial 
"trapped" second sound mode predicted by the theory leads to a sharp resonant dip in the reflected 
signal from an external second sound pulse propagated towards the interface when its horizontal 
phase speed matches that of the interface mode. The influence of the interface on thermal fluctua- 
tions in the bulk superfluid is shown to lead to slow power dependence of the order parameter, and 
other quantities, on distance from it. 

PACS numbers: 05.70.Jk, 05.70.Ln, 64.60.Ht, 67.40.Pm 



In a seminal paper, Onuki showed that when a uni- 
form heat current Q = Qz is driven through a sample 
of 4 He very close to the superfluid transition, a situation 
can occur where the induced temperature profile divides 
the system into separate upstream (z < 0) normal and 
downstream [z > 0) superfluid regions. Within the in- 
terface between them (taken as centered on z — 0), the 
primary mode of heat transport converts from thermal 
diffusion, with temperature gradient VT = — Q/ft where 
k is the thermal conductivity (diverging at the superfluid 
transition), to superfluid counterfiow, with a supercur- 
rent j s flowing toward the interface, and normal current 
j„ = — j s oc Q carrying the heat away from it. The latter 
flow shorts out all temperature gradients, leading to an 
asymptotically constant temperature Toa{Q) deep on the 
superfluid side, z — > oo @]. The interface width, £(Q), 
diverges as Q — > 0, and plays the role of the fundamental 
correlation length in the system. Correspondingly, the 
order parameter ip vanishes in the normal fluid, smoothly 
turns on through the interface, and takes the helical form 
ip(z) = \4>oo{Q)\ e ~ lk °°^ z deep in the superfluid, where 
the phase gradient fcoo — —mv s /h is proportional to the 
superfluid velocity. A mean field calculation of these pro- 
files is shown in Fig. ^ 

This initial work spawned a series of experimental |3( 
and theoretical Q investigations into this, and related, 
nonequilibrium superfluid critical phenomena Most 
relevant to the present work, in Ref. |(| the dynamics of 
the interface under external forcing was considered, and 
it was predicted that an interfacial second sound mode 
exists, in which perturbations travel along the interface 
as waves with a well defined sound speed c(Q), and higher 
order damping constant D(Q). The waves are trapped 
in the sense that their amplitude dies exponentially into 
the bulk phases on either side 6 J . In this paper an exper- 
iment is proposed, and the underlying theory developed, 
to verify the existence of this mode via second sound 
scattering. It is shown that when a pulse is propagated 
toward the interface from the superfluid phase (see Fig. 
EJ), strong resonant absorption occurs when its horizon- 
tal phase speed matches c(Q), leading to a sharp dip, 
with depth scaling with D(Q), in the reflected amplitude 
(Fig. |3 below) . In a related effect, it is shown that ther- 
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FIG. 1: Scaled mean field steady state mean field tem- 
perature, M, and order parameter magnitude, |*|, profiles 
through the interface as derived from Jjy, using parameters 
ao = bo = 1, do = 2 and co = 0, which produce Moo = —1.936, 
I = 1.216 and K^, = 0.676. 



mal excitation of these same modes leads to slow power 
law corrections in l/z to the order parameter and other 
quantities. 

The analysis is based on the Model F equations of 
Halperin and Hohenberg Q ■ These are derived from the 
Hamiltonian 



H = J d d r^-\V^ + \ 



r \^\ 2 +u \^\ 4 



;Xo L (m- xoh + 7oXoH 2 ) 2 



(1) 



where "0(r, t) is the superfluid order parameter, m(r, t) is 
a linear combination of mass and energy density, and the 
last term provides the crucial coupling between the two. 
The equations of motion are obtained from 
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where 9^ and 9 m are thermal white noise sources and 
W = Q[5(z — z\) — 5(z — Z2)] provides a source of heat 



at z\ — > — oo and a sink of heat at z 2 +00. The 
local temperature is defined as /i(r,i) = SH/Sm(r,t) — 
Xg 1 m + 7o|V'| 2 — ho- The various constant parameters 
may be partially determined by fits to experimental data 
The basic (mean field) length in the problem is Iq = 
(Ao/2xo7oQ) 1 ^ 3 , and it is convenient to define rescaled 
time and space variables R = r/Zo and r = t/to, with 
to = Z^/RelV The equations of motion take the form 
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where a = 5 /27oXo(Rer ), 6 = 7oXo5o/2A m , c = 
Imro/ReTo, d = 2u / X oli, * - ^o^e-^^, M = 
[ro + 27oXo/j]Z 2 , p- = lim^oo fi(z) is the asymptotic su- 
perfluid temperature, = limz-xx, M(Z), and J = 

Im(4'*VR4') is the supercurrent density. The rescaled 
thermal noise sources 0* = (2y / uoZo/Rer )e~ I£ ' t) ' i *#'!/, 
and &m = (27oXo^o/Ao)^m have covariances that diverge 
as l^~ d oc Q^f An expansion in 4— d is then required 

for a full theory in the small Q critical regime Q , but is 
difficult for the full interface problem. However, the ex- 
periment proposed here is less concerned with subtleties 
of noncquilibrium criticality, than with basic hydrody- 
namical properties of the interface that are more sharply 
defined and easier to investigate at larger Q. A corre- 
spondingly simpler theoretical approach will be taken in 
which Q is assumed large enough that the noise may be 
treated as a perturbation on the homogeneous equations 
□ 

Writing * = |*|e^ and U = [|*|, <j>, M], one may 
readily obtain a noise free steady state solution with an 
interface (see Fig.^J, denoted by U (Z), with M (Z) = 
-Z, Z^-oo, and Uo(Z) -» [jtfool.-ifoo^Mao], Z^ 
+00, with the constraints K 2 ^ + I'J'ool 2 + Moo — and 

W^oo|*oc| 2 = 1. 

Long wavelength, harmonic excitations of the inter- 
face are solutions of the form STJ = U(R, r) — Uo = 
(5U q ^(Z)e lcl ( X ' y '- 4 " T , and may be thought of as being 
driven by an external source on the plane Z — Z 2 with 
fixed frequency u> and transverse wavevector q. These 
satisfy linearized equations 



[to + Ud z + U(d 2 z - q 2 )]8\J^{Z), 

(4) 

with Z-dependent matrix coefficients which follow by 
straightforward linearization of @ about Uo, but whose 
exact form will not be required here. Deep on the super- 
fluid side, the solution is a sum of incoming and outgoing 
scattered waves, 



SV^(Z) = A- m [V-e 



iq z z 



«outV qu e 



(5) 



obtained from the asymptotic form of @ with the re- 
placement dz — > iq z - Here A m is the amplitude of the 
incoming excitation, and a ut(q, w) is the relative ampli- 
tude of the reflected wave. The wavevector components 
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FIG. 2: Scattering geometry: An incoming pulse with am- 
plitude Air, and spectrum narrowly centered on horizontal 
wavevector q and frequency uj approaches the interface at ve- 
locity vj. The reflected pulse, with relative amplitude a ou t, 
moves away at velocity vj~ , accompanied by an excitation 
of the interface itself with relative amplitude a at moving at 
speed C§ along the direction of q. When the resonance con- 
dition u)/q = Cg holds, a ou t — (see Fig. |3J and the pulse is 
strongly absorbed by the interface. 



g+(q, w) > q z (q_,w) are the two solutions to the second 
sound dispersion relation 

= (l + d )uj 2 +Aa K oo q z Lj- 2ag|* 00 |V (6) 
- 2ag(|^ 00 | 2 - 2Kl)q 2 z + 0(u 3 , ql, q z q 2 , . . .), 

where stability requires that |^oo| 2 > 2K Qa 
0. The corresponding eigenvectors are normal- 
ized to have unit (^-component, Vq W = [i(u> + 
2aoK 00 qf)/2a \V 00 \, 1, -tw/ao] with 0[lj 2 , (qf) 2 , q 2 , . . .} 
corrections in the first and last components. To the 
exhibited order, © may be put in the elliptical form 
{q z - q z ,c) 2 /(Aq z ) 2 + q 2 /(Aq) 2 = 1, with ellipse center, 
and semi-major axes given by q zfi — -Z^oow/ao(| v I'oo| 2 — 
2Kl), Aq z = q z ^(l + d )(\y oo f/2Kl-l) + l, 
Aq = Aq zy / 1 - 2Kl D /\^ 00 \ 2 . The direction of pulse 
propagation is determined by the group velocity 



v s(q^) 



1 + q z ,c(<lz - q z . c )/Aq 2 



(Aq) 2 ' (Aq z ) 2 



so that the v+ z > > v g z 



(7) 

indeed have opposite sign, and 
for each q, uj one may form pulses propagating towards 
and away from the interface (see Fig. [2J ■ 

The dissipative normal phase dynamics implies 
an exponentially decaying asymptote SMq^, = 

A in a N (q,u!)e iq * z , Z < 0, with Imgf < determined 
by the diffusion relation, i(bodo/2ao)oj = q 2 + {q^) 2 - 
In solving fll)!. we will formally treat u,q as 
small parameters of the same order, thus viewing 
Cs = ijj/q as a fixed 0(1) parameter. It follows 
that qf = qEZo«U 2k , = ? 1/2 Er=0^9 fc with 
known coefficients, 



-(1 



tt s ± Att V1-1/(AC)2, 



f i)^JbodoC s /^ao, 
where the ratios 



3 



7Tq = q z , c /q, Attq = Aq z /q, AC = Aq/q are functions 
only of Cs- 

The appearance of q 1 ! 2 in the normal phase leads one 
to expect the solution to Q to have an expansion in 
q 1 ' 2 rather than q or q 2 : SU^(Z) = J2Zo q k/2 SV k (Z), 
£*W,iV = Sfc!Lo o^ ut ' N q k ^ 2 . Asymptotic matching will 
produce simultaneous expansions in powers of q 1 ' 2 and 
Z (coming from the expansion of the exponentials) on 
each side of the interface, and matching each monomial 
coefficient will allow determination of the unknown coef- 
ficients a™ t,7V . Defining = Lo + Lidz + \j%d\, sub- 
stitution of this expansion into J3J leads to the sequence 
of relations 

L Z <SU = 0, L Z< 5U! = 

L Z <SU 2 = iC s L 3 SV , L Z SU 3 = tC s L 3 5V 1 

L Z 6\J 4 = iC s ~L 3 6U 2 -L 2 6U , (8) 

etc. Therefore, SXJq.i correspond to zero frequency per- 
turbations of the interface, i.e., infinitesimal motions that 
simply produce a new steady state. There are two of 
these: a global phase rotation v = [0,— /foo,0] (nor- 
malized using Kaa for later convenience), and a uniform 
translation of the interface vi = 8zUq. The asymptotes 
are v x v for Z » 0, 5M X -> -1 for Z < 0. The 
appearance of "Lz^k in (jSJ) at each order means that 
the solution can be determined only up to an arbitrary 
linear combination /ifcVo + v^wx whose coefficients must 
be determined from the matching conditions. 

To begin, one obtains <5Uo,i = /io.ivo + ^o,i v i- The 
O(q jf? 1 / 2 ) matching conditions yield v§ = —aff = 0, 
vi = —oti on the normal side, and 1 + af] 11 ' = — -R'ooMOj 
a° ut = — Xoo(mi + v i) on the superffuid side. 

At O(q) one obtains <5U 2 = 7/JoC s v 2 + + ^2^1, hi 
which V2 satisfies Lzv 2 = L3V0. It follows that v 2 is the 
change in shape of the interface under a perturbation of 
the heat current, Q — > 1 + SQ, and results from compres- 
sion and rarefaction of the heat current in the vicinity 
of the interface due to the incoming wave. Adjusting Q 
leads to a simple rescaling of Uo, and one obtains exact 
result v 2 = [K oo /2a M oo ){Zv 1 + [\V \,Q, 2M }) jfi], with 
asymptotes v 2 -> (K oo /2a M oo )[\^ oo \,-K oo Z,2M oo ], 
Z > 0, and M 2 -> -3K !>o Z/2a M oo , Z < 0. The most 
important matching condition now arises from the term 
linear in Z on the superfluid side, which produces the ad- 
ditional constraint a^Vo + = —fi CsK% /2aoM 00 , 
using which one obtains, 



-(1 + Koofia) = - 



tt _ K OD Cs/2a M o 



^ooCs/2a Moo' 



(9) 



along with v-i = —ot% , and v\ = 
3K 00 ^ioCs/2aoM 00 7rQ . The latter corresponds to 
the actual spatial amplitude of the interface motion, 
5Z(q,u) = A Ui v\q x l 2 . The numerator in © vanishes, 
implying full absorption of the incoming wave, for C's — 
C s ,o, with 



Cl ° ~ 2d Q M 2 



4a 2 ,|* 00 | 2 M£, 



|* 0O | 2 (2|* c 



K 2 



(10) 



corresponding precisely (the square of) the interfacial 
sound speed (with physical value c = Cs. oW T o) found 
in Ref. @. 

At order 0(q 3 / 2 ) one obtains SU 3 = iC s (/iiV 2 +i/iV3) + 
/-*3Vo + in which V3 satisfies L^V3 = L3V1. It fol- 

lows that V3 is the perturbation to the interface profile 
induced by a change SQ of the incoming heat current on 
the normal side only, with that exiting on the superfluid 
side unchanged: XJ(Z,t) — U (Z — vr) cx w 3 (Z — vt). 
This causes heat to build up behind the interface, moving 
it forward at an instantaneous speed v cx bQ)Z\. This 
heating effect leads to dissipation of the interface motion, 
which is singularly damped [6j , at rate cx q 3 ^ 2 rather than 
the bulk q 2 [which would arise from corrections to ©]. 

Although an analytic form for V3 is not available, it 
can be obtained numerically. For the purposes of deter- 
mining fii, i>2, a-2 1 ot° ut j only two numbers are required, 
namely the coefficients k 3 , m 3 in the asymptotic forms 
V3,$ — > k 3 Z, Z > 0, and v 3tM — * {bod / Aa )Z 2 + m 3 Z, 
Z <C (there may also be constant terms, but these 
may be absorbed into fi 3 ,v 3 , which remain undeter- 
mined at this order). In addition to the 0(q) constraint 
-.Koo(mi + "1)3 one obtains from the matching: 



a1 ut 



«/C s )a 2 v + (3if oo /2a M oo )/ii 
«/C s )«r + (^/2a M oo ) Ml 



v\m 3 
v\k 3 , 



(11) 



with solution 



3Cf (tt+ - 7r - ) k 3 + K 2 j2avM a 



2a M oo TT$ r 



..N 



-v 2 



Cs_ 



m 3 vi 



3« ut 



) C s /2a Mo O ) 2 



2a M a 



(12) 



Since 7Tq is complex, the zero of the first order combi- 



nation Q!out(?) ~ £ 
complex value to = Cs$q 
dissipation parameter 



+ g 1 / 2 a5 ut , is now shifted to the 
l-l)D Sfi q 3 / 2 + 0(q 2 ), with 



D s ,o 



9^ 0O (Cg) 7 / 2 (fc 3 + i^ 2 c /2a M oo ) 



(13) 



and corresponding physical value D — D^l^ 2 /tq. 

Using the parameters ao = bo = 1, do = 2, Co = 0, 
one obtains Cs = 1-089, k 3 = 0.488, L> s ,o = 0.143, and 
777,3 = 0.637 |6J. The resulting magnitude and phase of 
the reflection coefficient a^out), plotted as a function of 
Cs for various q, are shown in Fig. [3] 

Finally, consider the effects of thermal noise, under the 
assumption, as discussed above, that the scaled noise am- 
plitude is sufficiently small that it may be treated within 
the linear response regime. Using the mode decompo- 
sition U = J2q,u A ciA T ) SlJ ^(Z)e lCi < X ' Y \ one obtains 
an equation of motion for the amplitude 



[d T + icj)A qw (r) = 0q, w (r), 



(14) 



in which the white noise 0q >u ,(r) is the appropriate projec- 
tion of 0>j , 0m onto the mode eigenvector. The solution 
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FIG. 3: Magnitude and phase (inset) of the reflection coeffi- 
cient a ut = a§ ut + a ° ut \/9 plotted as functions of the hor- 
izontal phase speed Cs for q = 0,0.1,0.2,0.3,0.4,0.5, using 
scaled Model F parameters ao = bo = 1, do = 2, cq = 0. 
In both figures, sharper curves correspond to smaller q. The 
minimum allowed phase speed C™ m occurs for q^r — q z ,c, 
hence a purely horizontal group velocity. The large Cs asymp- 
tote of ao ut corresponds to a high frequency pulse propagated 
nearly straight at the interface. For small q the minimum oc- 
curs at Cs = Cs,o + Ds,oy/q with value B Ds,o^/q, B = 
2a 2 , Ml ( | * «, | 2 - 2 K 2 X ) /9 K 2 X j * «, | 2 (Cg ) 3 , scaling with the in- 
terfacial dissipation constant. The phase switches rapidly (in- 
stantaneously for q — > 0) between near-7r to near-zero as Cs 
increases through C s . 



(15) 



(where the small dissipative negative imaginary part that 
gets added to cu at order q 2 is actually required here to 
ensure convergence) inserted back into U, allows one to 
compute stochastic averages of various quantities (cal- 
culations are tedious and badly encumbered by matrix 
indices, and will be presented in detail elsewhere). Not 



surprisingly it is the phase correlations that are the most 
important, decaying at long distances as a power law 
(<^(R)<50(R')) oc |R - R'| 2 " d (with d = 3 here) deep 
in the superfluid phase, with a complicated anisotropic 
coefficient depending on the angle between R R' and 
the heat flow direction z. The extra 0(q) factors in the 
|\&|,M components of Vq U produce weaker power laws 



IR-R'I 



also 



(<J|*(R)|J|tf (R')|), {SM{K)SM(K')) cx 
with anisotropic coefficients. 

Closer to the interface, there are residual correlations, 
coming from interference between the incoming and re- 
flected waves in V^,, that produce Z-dependence in 
the local fluctuations: (Scj)(oo) 2 ) - (S4>(R) 2 ) cx l/Z d - 2 , 
(SM(oo) 2 ) - (SM(R) 2 ), (<5|*(oo)| 2 - <5|*(R)| 2 ) cx l/Z d , 
where the argument oo is shorthand for Z — > oo. The 
phase fluctuations have a very strong effect on the com- 
plex order parameter: 



(*(R)) 



(|*(R)|)(e^ (R) ) 
*o(R)e"^ <l5<!i(oo)2 



: * (R)e-^^( R ) 2 > 
(l-BZ 2 - d ), (16) 



where B is a coefficient, and the Gaussian property of the 
noise has been used to average the exponential. There 
are two interesting effects here: first, phase fluctuations 
can significantly reduce the magnitude of the order pa- 
rameter even in the linear response regime [validity of 
(fTfi)) requires only the weaker assumption that the space- 
time derivatives of 8<p, as opposed to 6<f> itself, be small]. 
Second, the presence of the interface induces a slow 
power law approach of |(W(R))| to its asymptotic super- 
fluid value, contrasting with the exponential approach of 
'I'o(R)- This power law is not induced by the positional 
fluctuations of the interface itself, which remain strongly 
bounded 0, but by the effects of its mere presence as 
a reflecting boundary on the long-range bulk superfluid 
correlations. 

At linear order, the average temperature (M(R)) = 
Mq(R) remains equal to its mean field value, with the 
power law visible only in the variance. However, it is 
likely that similar power laws will be induced in (A/(R)) 
if nonlinear corrections are taken into account. 
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